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Connecting the dots: from theory to (gravitational wave) observations 


Connecting 
theory and 
observations 


Credit: H. O. Silva 


* Waveform templates crucial for signal characterisation 

* Formation and (stability) properties of black holes and compact objects 

* Construct consistent inspiral-merger-ringdown waveform models in and beyond General Relativity 
* Combine theory-agnostic waveforms (e.g. Maggio, Silva, Buonanno, Ghosh ‘22) and theory-specific waveforms 


See Laura and Mohammed’s talks for PN calculation; Geraint's and Nathan's talks on parameter estimation. 


Landscape of Numerical Relativity beyond GR 


In four spacetime dimensions the only divergence-free symmetric rank-2 tensor constructed 
solely from the metric and its derivatives up to second order, and preserving diffeomorphism 
invariance, Is the Einstein tensor plus a —— constant. [Lovelock '69-'72] 
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Focus today: higher derivative gravity 


* Higher curvature corrections 
relevant in strong-curvature regime 
e Bottom-up: series expansion in curvature 
e Top-down: inspired by low-energy limit of quantum 


gravity candidates, e.g., string theories 
(Gross & Sloan '87, Kanti et al '95, Moura & Schiappa 06, Cano et al '21) 
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Scalar Gauss-Bonnet Dynamical Chern-Simons Effective field theory of gravity 


G = Rapea R” — 4R pR” + R* 


How to formulate higher derivative gravity theories? 


* Field equations may have (i) mixed hyperbolic-elliptic character, (ii) higher time derivatives,... 
* Causal and well-posed initial value formulation necessary for numerical stability 


Modified gauge conditions 


Fixing a la Israel-Stewart 


Order-by-order expansion 
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See Aron's and Miguel's talks 


Building a black hole 


Building 
a (hairy) black hole 


Ultralight scalars: E.g.: in quadratic gravity: 
Black holes anchor long-lived scalar clouds Black holes have scalar charge, but neutron stars do not! 
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See Thomas’ talk for solution space, Nico’ talk for QNMs + mode stability 


Growing monopole and dipole hair in quadratic gravity 
(A. Hegade, E.R. Most, J. Noronha, HW, N. Yunes '21, '22) 


How does the black hole hair form? 


* simulated stellar collapse 


* scalar field growth triggered 
by formation of event horizon 
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Scalar clouds non minimally coupled to gravity 
(C. Richards, A. Dima, HW ‘23) 


* Evolving massive scalar cloud around single black hole massive dCS hair 
"WE" r i . . scalar cloud 
e Non minimally coupled to gravity via Pontryagin density 


c 9+- S RR - 420 2 0 el + 


(2) 
B ele 


Massive scalar cloud Oscillating Chern-Simons hair 
t-2500M 46, 0.. 
0.006 20 0.006 
0.004 E 0.004 
10; 
0.002 0.002 
D 
x 0000 = o 0.000 
> N 
—0.002 —0.002 
ae 15| a -0.004 BH Spin: a/M = 0.99 
-0.0006  -20|uM-0.42 P Particle mass 


—20 —10 Mu — 0.42 


X/M 


See: massive dCS scalar around Schwarzschild (Macedo '18), Superadiance of massive scalar around dCS black hole (Alexander et al '22) 


Binary black hole mergers 


Merging 
black holes beyond GR 
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Canuda 


Open-source numerical relativity code 
for fundamental physics 


(9 Witek, Helvi; ®© Zilhao, Miguel; @ Bozzola, Gabriele; @ Cheng, Cheng-Hsin; ($ Dima, Alexandru; @ Elley, Matthew; © 
Ficarra, Giuseppe; ($ Ikeda, Taishi; @ Luna, Raimon; @ Richards, Chloe; @ Sanchis-Gual, Nicolas; (9 Silva, Hector 


https://bitbucket.org/canuda/ 
© Canuda Numerical Relativity 


Coalescence of black holes with scalar hair ( charge”) 


e scalar dipole radiation 
e Gravitational wave phase shift 
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PN Waveform (Shiralilou, Hinderer, Nissanke, Ortiz, HW '20, '21) 
[See also: Yagi+ '11, Julie+ 18, '19, Lyu+ '22] 


[Equivalent numerical approach in dCS: Okounkova et al '17-'19| 
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Scalar Gauss-Bonnet gravity Il: dynamical (de-)scalarization 
(Silva, Elley, HW, Yunes '21, '22) 


* Consider quadratic coupling of scalar field (0O -— u? \ = 0 u? ~ le 
— ! " eff eff 
e Scalarization of single black holes if ps < 0 


* Coalescence “""® dynamical scalarization or desclarization 
e Nonlinear phenomenon may affect only merger and postmerger! 


time = 492 M (iteration = 100096) 


[Single BHs: Silva et al '17, Doneva et al '17, Antoniou et al '17, Blazquez-Salcedo '18, Macedo et al '19, Ripley & Pretorius '20, 
opin-induced scalarization: Dima et al '20, Hod '20, Doneva et al '20, Herdeiro et al '20, Berti et al '20, 
Binaries: Silva, HW, Elley, Yunes '21; East, Ripley '21; Review: Doneva et al ’22] 11 


Scalar Gauss-Bonnet gravity w/ modified generalized harmonic gauge 
(Corman, Ripley, East '22) 
l j AGB 
e Scalar Gauss-Bonnet gravity Z = R — Pi V o) v - 


* Modified generalized harmonic gauge 
* Breakdown near merger for some parameters 


Mass ratio g = 1 Mass ratio g = 1/2 
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[See also Kovacs & Reall '20, East & Ripley '21, '22] 12 


Four derivative scalar-tensor theory 
(Sole, Clough, Figueras '22) 


* General scalar-tensor theory w/ up to four derivatives 
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Fourth order effective field theory of gravity 


(Cayuso, Figueras, Franca, Lehner '23) 


* Effective field theory of gravity (parity even and in vacuum) 
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L = R— E — 
* Using 

i. order reduction 

ii. fixing equations a la Israel-Stewart 
e Binary black holes w/ q=1, d = 12M 
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|J. Cayuso, Ortiz, Lehner '17, Allwright & Lehner '18, Franchini et al '22, Lara et al '22, Gerhardinger et al ’22] 
Note: order reduction refers to inserting the eoms for the Ricci tensor, thus shifting its contribution to higher-order. 
Sometimes this terminology is used for the order-by-order expansion, but not here 


A note of caution 


What about degeneracies with 
e Unmodeled GR effects? 
e Environment? 


Black holes binaries in dark matter environments 
(C.-H. Cheng, Ficarra, HW in prep.) 


* Consistent evolution of massive scalar cloud around binary black hole | j— u^ o = 0 
* Goal 1: new gravitational wave probe for new particles 
e Goal 2: toy model for environmental effects during merger 

e Constraint satisfying initial data and eccentricity reduction 
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[See also Bamber et al '22] 5 / M t —0M 


Summary 


Status update in Numerical Relativity beyond GR (focus: higher derivative gravity) 


* Proof-of-principle simulations of compact binary coalescence in classes of beyond GR model 
e Scalar charges typical yields scalar (dipole/quadrupole) radiation 

e Dephasing of gravitational wave signal 

e Surprisingly simply merger signal — not as nonlinear as may have been expected 


Where do we go from here? 
e Synergy w/ PN and constructing full IMR models -> Waveform catalogs? 
* Dictionary with theory-agnostic models for the plunge+merger (e.g. EoB) 
* Concerns? 
* Hobustness: waveform comparison and benchmark from different approaches and codes 
* Eccentricity reduction of initial data 
* Degeneracies with unmodelled GR or environmental effects 


Thank you! 
E ILLINOIS 
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